arXiv:1503.06624vl [quant-ph] 23 Mar 2015 
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The quantum Fourier transform (QFT) is a key ingredient of several quantum algorithms and a 
qudit-specific implementation of the QFT is hence an important step toward the realization of qudit- 
based quantum computers. This work develops a circuit decomposition of the QFT for hybrid qudits 
based on generalized Hadamard and generalized controlled-phase gates, which can be implemented 
using selective rotations in NMR. We experimentally implement the hybrid qudit QFT on an NMR 
quantum emulator, which uses four qubits to emulate a single qutrit coupled to two qubits. 


I. INTRODUCTION 

The size of a quantum register can be increased ei¬ 
ther by increasing the number of qubits (in the standard 
model of a quantum computer) or by increasing the num¬ 
ber of accessible logical states in each quantum element 
i.e by using qudits (quantum digits) [iH^- For qudit com¬ 
puting to become a reality, it is essential to concretize 
the theoretical framework for computation in terms of a 
set of quantum gates that are universal and can imple¬ 
ment any algorithm efficiently. It is however not always 
straightforward to generalize a qubit gate to the qudit 
scenario and in some cases, more than one kind of gen¬ 
eralization is possible, each retaining some features of 
the original qubit gate @. A gate library for qudits has 
been proposed, containing one-qudit and two-qudit gates 
such as the CINC gate, generalized SWAP and general 
controlled-X (GCX) gates Several qudit comput¬ 

ing proposals using different physical systems have been 
designed to incorporate “hybrid” qudits of different di¬ 
mensions [iMl , wherein multi-qudit gates are general¬ 
ized to transform two or more qudits of different dimen¬ 
sions. A set of hybrid quantum gates have been designed 
to act on qudits of different dimensions, including the hy¬ 
brid SUM, SWAP, Toffoli and Fredkin gates [ll,[i^. Qu¬ 
dits have been used as a resource in secure quantum com¬ 
munication llTII , quantum key distribution [l8| , quantum 
computing [l9| | and quantum cloning [^ . 

NMR has been a fruitful testbed for implementing 
quantum information processing ideas and is now increas¬ 
ingly being used to manipulate qudits. NMR qudit imple¬ 
mentations include using quadrupolar nuclei oriented in 
a liquid crystal for information processing, studying 
the dynamics of nonclassical correlations and finding 
the parity of a permutation using a single qutrit [2^ and 
a single ququart [^. Logic gates and pseudopure states 
for a ququart have been implemented using ^^Na and ^Li 
nuclei oriented in a liquid-crystalline matrix psi . . 
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The quantum Fourier transform (QFT) needs to be 
specifically tailored for qudit computers as it plays a 
key role in several quantum algorithms such as factor¬ 
ization [ 23 , quantum phase estimation [ 2 ^ and the hid¬ 
den subgroup problem [ 2 ^. An efficient decomposition 
of the QFT for qubits achieves an exponential speedup 
over the classical fast Fourier transform [s^, and has 
been experimentally implemented by several groups on 
different physical systems [3ll - l35l| . Decompositions of the 
QFT fo r a udit systems have been worked out by several 
groups 

In NMR systems with “strong” spin-spin coupling 
terms in the Hamiltonian, the identification of spins as 
qubits becomes problematic. In such cases, quantum 
gates have been implemented by selectively manipulat¬ 
ing individual transitions between pairs of energy lev¬ 
els [iO, |4l|. Such “effective” qubits have been con¬ 
structed by a mapping to the non-degenerate energy 
eigenstates of a single qudit and have been experimen¬ 
tally demonstrated by partially orienting a ^^^Cs nuclear 
spin (S=7/2) in a liquid-crystalline medium (43 - l43| . We 
carry these ideas forward by carving out a system of two 
qubits and a qutrit (henceforth referred to as a QQT 
system) out of four qubits (henceforth referred to as a 
QQQQ system). Specifically we use the four coupled 
qubits of 5-Fluorotryptophan to emulate a hybrid system 
of two qubits and a qutrit (a 2 (8)2 (8) 3 dimensional Hilbert 
space). The qubit architecture of 5-Fluorotryptophan 
due to its specific spin-spin coupling strengths and molec¬ 
ular symmetries, gives rise to twenty lines in the NMR 
spectrum instead of thirty two lines expected for the non¬ 
overlapping energy levels of four qubits. This matches 
the spectrum expected for a system of two qubits cou¬ 
pled to a qutrit. This system obviously cannot be used 
as a full-fledged QQQQ system of four qubits due to the 
degeneracy of transitions. However, by re-labeling and 
creating a mapping of corresponding energy levels, we 
are able to show that this QQQQ system can indeed em¬ 
ulate a QQT system and that the allowed non-degenerate 
transitions are sufficient to implement the required quan¬ 
tum gates. 

In this paper we describe a decomposition of the hy¬ 
brid qudit QFT, using a set of generalized Hadamard 
and generalized hybrid controlled-rotation gates. The 
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scheme is general and can be implemented on any physi¬ 
cal hardware for a hybrid qudit quantum computer. We 
implement the hybrid QFT on a QQT system (using the 
four qubits of 5-Fluorotryptophan as a QQT emulator). 
The requisite hybrid qudit gates can be realized in NMR 
using selective rotations. Partial state tomography of 
the final state was achieved by a set of 19 specially de¬ 
signed experiments to reconstruct desired portions of the 
Hilbert space. The results of the tomography demon¬ 
strate the success of the mapping of the QQT system 
onto the QQQQ system as well as the implementation of 
the hybrid qubit-qutrit QFT on the system. 

The material in this paper is organized as follows: Sec¬ 
tion mi describes the basic gates required for quantum 
computing with qudits and the circuit decomposition of 
the qudit QFT. Section mil describes the experimental 
implementation of a hybrid qubit-qutrit QFT on a four- 
qubit NMR quantum emulator. Section lTVl contains some 
concluding remarks. 


controlled-rotation gate R^f. (with j being the control 
qudit and k the target qudit) is given by P,[l3: 


■ ■ ■ Xj ■ ■ ■ Xk - ■ ■ XN-i) = 


exp 


— 2'KiXjXk 
. Ylp=j'^p , 


\xo - ■ ■ Xj ■ ■ ■ Xk - ■ -XN-l) 


( 4 ) 


where dp is the dimension of the qudit; k > j has 
been assumed here, in order to explicitly define the ac¬ 
tion of the gate on the N-qudit basis state, however inter¬ 
changing j and k does not alter the gate operation. For 
dp = d = 2, the hybrid two-qudit controlled-rotation gate 
reduces to the standard two-qubit controlled-rotation 
gate Rj^k- 


B. Hybrid Qudit QFT 


II. QFT DECOMPOSITIONS 


The action of the QFT on the basis vectors of a hybrid 
qudit system is given by [T^ Is^ : 


The state of a system of N hybrid qudits each of a 
different dimension dp {p = I ■■■ N) can be written in 
terms of an orthonormal basis of product states 

Is) = |so) 0 |si) 0 • • • 0 |sAr_i) 

= \X0---Xj---Xk---XN-l) ( 1 ) 

where xj G {0,1,... dp}. For dp = d = 2, this reduces to 
an N qubit state, with x = ^'^*4 Xj G {0,1} 

being binary integers. However, x does not retain this 
simple form for hybrid qudits with dp > 2. 


A. Qudit Gates 


The Fourier gate (Fp) is a single-qudit gate that creates 
a superposition of all basis states of the qudit, with equal 
amplitudes, with its action on the qudit in an 7V-qudit 
system given by P,|3^: 


Fp\xj) = 

[0,1,2,... (d-1)] (2) 

For d = 2 (qubit) and d = 3 (qutrit) the Fourier gate 
reduces to the Hadamard gate (H) and the Chrestenson 
gate (C) respectively: 
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For a system of N “hybrid” qudits, each of different di¬ 
mensions di, (i 2 , ■ ■ ■ dN, the action of a two-qudit hybrid 



where the states \y) have the same form as jx). 

The circuit for implementing the QFT on N hy¬ 
brid qudits can be decomposed as a set of single-qudit 
Fourier gates interspersed with two-qudit controlled- 
rotation gates: 

PiRi,2Ri,3 - ■ ■ Ri,nF2R^,3 - ■ ■ Fn-iR^_i^nFn. For N 
qudits each of a different dimension d^, the total dimen¬ 
sion of the Hilbert space is D = Jl^i ^'^‘4 it has to 
be kept in mind while implementing hybrid controlled- 
rotation gates that the dimensions of the control and 
target qudits may be different. After the implementa¬ 
tion of the QFT using the above decomposition, the bit 
values of the resultant state appear in the reverse order. 
A sequence reversal can be achieved either by applying a 
series of multi-valued permutation gates @ or by reading 
off the result in the reverse order. 


III. NMR IMPLEMENTATION 

A. The hybrid QQT system 

We consider a QQT system with a Hamiltonian gov¬ 
erned by m, 

= Wi/i, + W2/2. + + DgiSll - / 2 ) 

+ Jl2llzl2z + J23I2ZI3Z + Ji^hzhz ( 6 ) 

where 1, 2 and 3 label the qubits and the qutrit respec¬ 
tively, liz is the z-component of the magnetization vector 
of the qubit (qutrit), P = * = x,y,z is the 

total magnetization of the qutrit, and a;i,a; 2 ,H 3 denote 
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FIG. 1. Energy level diagrams: (a) For a QQQQ system 
showing 16 energy levels corresponding to four coupled spin- 
1/2 particles. Spin states are given below each energy level, 
0 and 1 correspond to spin states and — i respectively. 
There are 32 transitions originating from single spin flips. 
Transitions belonging to spins 1 and 2 are shown in bine and 
brown respectively while transitions beloirging to spins 3 and 
4 are shown in purple, (b) For a hybrid QQT system showing 
12 energy levels and 20 transitions with energy levels labeled 
in the computational basis; 0 , 1 and 2 correspond to the 
qutrit spin states -1-1, 0 and —1 respectively. Each qubit un¬ 
dergoes six single spin flips shown in blue (spin 1) and brown 
(spin 2) respectively. Qutrit transitions are shown in purple. 


their Larmor frequencies i.e. the free precession of the 
qubits(qutrit) in a static magnetic fieid. Dq represents 
the effective quadrupoiar spiitting of the qutrit ieveis and 
the corresponding term in the Hamiitonian accounts for 
the static first-order quadrupoiar interaction. JijS denote 
the strength of the scaiar coupiing interactions between 
and qubits, and Jb represents an interaction term 
invoiving the qutrit. 

The QQT Hamiitonian in Eqn. (jH]) can be emuiated by 
four qubits such that two of the qubits mimic a three-ievei 
system. We use spins 3 and 4 of our four-qubit system to 
mimic a qutrit, such that the qutrit subspace containing 
ieveis 1,2 is emulated by the third qubit and the qutrit 
subspace with levels 2,3 is emulated by the fourth qubit, 
with the corresponding chemical shifts: + Dq 

and W4 = Q 3 — Dq. 

The Hamiltonian of a four-qubit QQQQ system (as¬ 


suming a specific coupling pattern J 13 = J 14 = Jj^g, 
J 23 = J 24 = J 23 J 34 = 0 ) is given by 

ffQQQQ — loiliz -f UJ 2 I 2 Z + ^ 3/32 -f a;4/42 + Ji2lizl2z 

+■^ 23 ^ 22(^32 + liz) + Jishzilsz + liz) (7) 

where liz represents the z-component of the magnetiza¬ 
tion vector of the spin. The NMR spectrum and mul- 
tiplet pattern of ijQQQQ with the assumed J coupling 
values resembles the first order NMR spectrum of a QQT 
system (wherein each qubit multiplet contains six transi¬ 
tions and the qutrit multiplet contains eight transitions). 
Figure 1 shows the energy level diagram of the four-qubit 
QQT emulator (energy levels numbered from 1-16) and 
an actual QQT system (energy levels numbered from 1- 
12 ). 

We use a 5-Fluorotryptophan molecule to emulate the 
QQT system, with one and three ^H nuclei repre¬ 
senting the four qubits. The molecular structure and 
NMR parameters are depicted in Figure 2(a) and the 
thermal equilibrium ID NMR spectrum is shown in Fig¬ 
ure 2(b). The Ti and T 2 relaxation times for all four 
qubits range between 1.17 - 3.96 s and 0.82 - 2.28 s re¬ 
spectively. The line intensities in the coupling pattern 
of 5-Fluorotryptophan differs slightly from that of the 
ideal QQT emulator, however hybrid qudit gates imple¬ 
mentation can be emulated on this system by tailoring 
the J—evolution intervals to obtain the desired angles 
of rotation. The multiplet pattern and single quantum 
NMR transitions of the QQT emulator are shown in Fig¬ 
ure 3; labels below each spectral line refer to the corre¬ 
sponding transitions as per Figure 1(a). The mapping 
between the transitions (labeled Ql, Q2, T) of the QQT 
system and the four-qubit QQT emulator (labeled ql, 
q2, q3, q4) is given in Table 1. For an ideal QQT emula¬ 
tor, the lines of the qubit occurring at uJi ± Jis =F Ju 
overlap with each other, while in the 5-Fluorotryptophan 
molecule, Jis yf Ji 4 {i € [ 1 , 2 ]) and the lines occurring at 
u)i ± Ji 3 =F Ji 4 are separated by 2 \ Ja — Ja j. However 
this molecule can emulate a QQT system whenever the 
intensities of the transitions of the spin occurring at 
uii ± Ji 3 =F Ji 4 are equal. 


B. Experimental implementation of the QFT 

The NMR pulse sequence for the experimental imple¬ 
mentation of a hybrid qudit QFT on the four-qubit QQT 
NMR emulator is shown in Figure 4. A hybrid controlled- 
rotation gate (where ‘j’ is the control-qubit and ‘fc’ 
is the target qutrit) can be realized by a z-rotation of the 
qutrit state by an angle (j S [0,1], k G [0,1, 2]). 

1 lm=j 

For the two qubits mimicking the qutrit, this controlled- 
rotation can be obtained by rotating both the qubits by 
the same angle with respect to the control qubit. Chemi¬ 
cal shift offsets are chosen such that spins 1 and 2 always 
remain on resonance. 
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TABLE I. Mapping between the transition frequencies of a 
QQT system (labeled Ql, Q2 and T) and a fonr-qnbit QQQQ 
system (labeled ql, q2, q3, q4) which acts as QQT emulator. 


Ql 

ql 

Q2 

q2 

T 

q3, q4 

1-7 

1-9 

1-4 

1-5 

1-2 

1-3, 2-4 

2-8 

2-10 & 

3-11 

2-5 

2-6 & 

3-7 

2-3 

1-2, 3-4 

3-9 

4-12 

3-6 

4-8 

4-5 

5-7, 6-8 

4-10 

5-13 

7-10 

9-13 

5-6 

5-6,7-8 

5-11 

6- 14 & 

7- 15 

8-11 

10- 14 & 

11- 15 

7-8 

9-11, 10-12 

6-12 

8-16 

9-12 

12-16 

8-9 

9-10, 11-12 

- 

- 

- 

- 

10-11 

13-15, 14-16 

— 

— 

- 

- 

11-12 

13-14, 15-16 


The action of the QFT in extracting state periodicity 
can be equally well demonstrated on an initial thermal 
equilibrium state as on a pseudopure state and we hence 
implement the QFT on a thermal equilibrium state [Slj. 
Beginning with the thermal equilibrium density operator, 
the first (f) pulse implements a Hadamard gate on the 
first qubit. The axis of rotation rii = cos^ii + sin^iy 
(aligned at an angle of 9i = ^ ^ with the x-axis) was 

chosen to obtain a 2 —rotation of the first qubit by an 
angle thereby circumventing the high time cost in¬ 
volved in implementing 2 —rotations. At this point in the 
pulse sequence, the qubits 2, 3 and 4 do not have any co¬ 
herence, and thus single spin 2 —rotations are redundant 
in their case. The first evolution period incorporates 
the controlled-rotation gates i?i ^2 and for a time 

Tij = (6ij is the desired rotation angle). The evo¬ 
lutions under the scalar coupling J 12 by an angle un¬ 
der the scalar coupling J 13 by an angle and under the 
scalar coupling J 14 by an angle ^ are achieved simultane¬ 
ously during the first evolution period. Refocusing pulses 
are introduced at appropriate points to obtain the desired 
rotation angles. The next module of the hybrid QFT be¬ 
gins with a (§)j^ pulse that implements a Hadamard 
gate along with producing an effective 2 —rotation of ^ 
on the second qubit (the axis n 2 is chosen to make an 
angle of 6*2 = f + ^ with the x-axis). This is followed by 
an evolution under the couplings J23 and J24 during the 
second evolution interval. The effect of evolution on the 
first spin is refocused by a pulse in the middle of the evo¬ 
lution period. The final Hadamard gate on the third and 



4129.m2 J 12 = 9.5H2 

-70622.2Hz J 13 = 2.5Hz 

4345 . 2 H 2 Ji 4 = 9.1Hz 

4381.2Hz J 23 = 9.9Hz 

J 24 = 4.6Hz 
J 34 = OHz 


(b) 

Spin 4 Spin 3 Spin 1 Spin 2 



■I.I.I.I . . . . 'I.I , . . , ' ' I ' ' I ' ' 

7.32 7.29 7.26 7.23 6.90 6.87 -125.01 


Wff (ppm) ujp (ppm) 


FIG. 2. (a) Molecular structure of 5-Fluorotryptoplian with 
qubits of interest shown in the rectangular block and num¬ 
bered from 1-4. Three nuclei and one ^®F nucleus form the 
four-qubit system. The chemical shift values (in ppm) and in¬ 
direct coupling constants (in Hz) are tabulated alongside, (b) 
NMR spectrum of the thermal state of 5-Fluorotryptophan. 
The and ^®F spectra are shown on the same scale. Colors 
of the spectral lines correspond to the spin color given in the 
molecular structure. The spins resonate at the frequencies 
uji, U 12 and 0 ) 3 , with Jij corresponding to the scalar coupling 
constant between spins i,j. 


fourth qubits is performed by a pulse (which ef¬ 

fectively mimics a single-qutrit Fourier gate). All pulses 
on qubit 1 are Gaussian pulses of 13 ms duration and 
the simultaneous excitation of qubits 3 and 4 is achieved 
using a 7 ms shaped pulse {Seduce. 100). All the pulses 
on qubit 2 are hard pulses. Since qubits 3 and 4 are only 
36 Hz apart it is difficult to individually address them in 
an experiment, and hence the controlled-rotation gates 
on these two qubits have been implemented without ap¬ 
plying individual refocusing pulses on them. 


C. Tomographic reconstruction 

Partial tomography of the QQQQ system was per¬ 
formed in order to reconstruct the desired QQT density 
matrix. First, a set of operators were designed to recon¬ 
struct the full QQT density operator in dxd dimensional 
Liouville space (in this case d = 2 ® 2 0 3). This set of 
operators was then extended to the analogous set of to¬ 
mography operations in a 16 x 16 dimensional operator 
space of the four-qubit NMR emulator. 

The complete characterization of a 12 x 12 dimensional 
QQT density operator requires the determination of 143 
variables. The 11 diagonal elements (populations) are 
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Qubit 1 (ixJi) Qubit 2 ( 102 ) 



(^) Qubit 3 (ixJs) 


/\ 

/ \ 



<n fo Tf< 


C*^) Qubit 4 (uoi) 




T12 


^"24 


FIG. 4. Pulse sequence for a hybrid QFT on a qubit-qubit- 
qutrit (QQT) system implemented on a four-qubit emulator. 
Flip angles and axes of rotation are displayed over each pulse. 
Rotation axes ni and n 2 are aligned at angles 9\ = f+ ^ and 
02 = f ^ with respect to the x-axis. Thick black rectangles 
represent selective pulses while thin black rectangles represent 
non-selective pulses; njS represent time intervals of evolution 
under the scalar couplings Jij. 


FIG. 3. Multiple! pattern: (a) Of qubit 1, (b) Of qubit 2, 
(c) Of qubit 3 and (d) Of qubit 4 of the four-qubit NMR em¬ 
ulator. Each multiple! is centered around its chemical shift 
value (black lines representing cji, ai 2 , ws and loa respectively), 
which is further split by scalar coupling interaction. The ratio 
of coupling constants J 12 : J 13 : ^ 14 , J 12 : J 23 ’■ J 24 , J 13 ■ J 23 
and Ji 4 : J 24 is the same as in the 5-Fluorotryptophan 
molecule. The line thickness is proportional to the spectral 
line intensity. Labels below each line refer to the correspond¬ 
ing transition as marked in Figure 1(a). 


obtained by applying an appropriate z-gradient to kill off- 
diagonal elements, followed by spin-selective rotations to 
project the diagonal elements onto experimentally mea¬ 
surable parts of the density matrix. The remaining 132 
elements are obtained by a set of 19 operators: III', YII', 
XIT, lYT, IXI', IIY', lYY', IXY', YIY', XIX', XXI', 
YYT, XYI', YXT, XYY', XXY', YXX', YYX' and IIA'- 
Grad^-IIT', where I is the identity operator, X(Y) refers 
to a single spin operator and primed operators corre¬ 
spond to qutrit operations. These operators can be im¬ 
plemented by applying the corresponding spin-selective 
^ pulses (or a no-operation for the identity operator). 
The last tomography experiment consists of a z gradient 
pulse sandwiched between two y pulses represented by A 
and T applied on the qutrit, of flip angles (—f) and {j) 
respectively. 

A corresponding set of 19 operations can be used to 
determine the off-diagonal elements of the four-qubit 
QQQQ system: IIII, YIII, XIII, lYII, IXII, IIYY, lYYY, 
IXYY, YIYY, XIXX, XXII, YYII, XYII, YXII, XYYY, 
XXYY, YXXX, YYXX and IIAA-Grad^-IITT. The last 
tomography experiment on the QQQQ system consists of 
a z gradient pulse sandwiched between two y pulses rep¬ 
resented by A and T, applied simultaneously on qubits 3 
and 4 of flip angles (—f) and (|) respectively. 


The theoretically expected and experimentally ob¬ 
tained tomographs of the real and imaginary parts of the 
density matrices obtained after applying a hybrid QFT 
on the four-qubit QQT NMR emulator are given in Fig¬ 
ure 5. On visual inspection, one can see that the general 
pattern of the weights in the density matrices is quite sim¬ 
ilar, indicating the success of the QFT implementation. 
Quantitatively, the fidelity of the reconstructed state was 
computed using the measure [3l|: 

p _ _ ^^(PtheoryPexpt) _ 

^(^^(pLoryPth eory ) ) \/(^^ ( Pexpt Pexpt) ) 

where ptheory and Pexpt denote the theoretical and ex¬ 
perimental density matrices respectively. The fidelity 
of the experimentally reconstructed density matrix was 
computed to be 0.83. The loss in fidelity of the exper¬ 
imentally reconstructed state can be attributed to sev¬ 
eral experimental factors, namely, rf pulse imperfections, 
imperfect refocusing of chemical shift evolution of un¬ 
wanted coherences during the evolution intervals, rf field 
inhomogeneity, and deleterious effects due to spin relax¬ 
ation. Our experiments point the way to an interesting 
direction of research in quantum computing with NMR 
namely, that of using complex coupled qubit topologies 
to emulate hybrid and higher-dimensional quantum sys¬ 
tems. 


IV. CONCLUDING REMARKS 

It has been proposed that qudit-based quantum com¬ 
puters are able to better optimize Hilbert-space dimen¬ 
sionality and are hence expected to be more powerful 
than the standard models of qubit quantum computers. 
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Re Im 



Re Pexpt Im 


One of the key quantum circuits in several quantum al¬ 
gorithms is the QFT and it is hence important to look 
for decompositions of the QFT specifically designed for 
qudits and hybrid qubit-qudit systems. We used four 
NMR qubits to emulate a hybrid system of two qubits 
coupled to a qutrit and implemented the QFT on this 
hybrid quantum computer. It is expected that these ex¬ 
periments will pave the way for the implementation of 
full-fledged qudit-based quantum computing proposals. 




FIG. 5. Tomographs of the real and imaginary parts of the 
(a) theoretically expected (ptheory) and (b) experimentally ob¬ 
tained (pexpt) final density operator obtained after applying 
a hybrid QFT on the thermal equilibrium state of the four- 
qubit QQQQ system used as a QQT emulator. The rows and 
columns encode the computational basis in binary order, from 
|000') to |112'). 
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